Incoherent tunneling surface diffusion for interacting adsorbates 
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The incoherent tunneling surface diffusion of light atoms is analyzed in terms of a Chudley-Elliot 
master equation constrained to nearest neighbors for moderate to high temperatures (above the 
crossover surface temperature). Recent helium spin-echo measurements for diffusing H and D atoms 
on a Pt(lll) surface are thus examined within the framework of the so-called bounce technique. 
A good agreement is found between theory and experiment, with the substrate friction coefffcients 
obtained being of the order of those reported in the literature for this kind of systems. Furthermore, 
following this model, an increase of the tunneling rate is predicted for decreasing surface coverages. 
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Hydrogen diffusion on metal surfaces has being the 
subject of an intensive applied and fundamental research 
for many years [ll-Q- The involvement of this process 
in the design of fuel cells or in the storage of H2 makes 
it of technological interest. Therefore, at a fundamental 
level, a good understanding of this dynamics is essential, 
for it provides important information on the adiabatic 
interaction between adsorbate and substrate. Indeed, 
because atomic hydrogen is the simplest element that 
may undergo chemisorption, its motion on a weakly cor- 
rugated surface constitutes a benchmark to study two 
possible competing diffusion processes: classical (acti- 
vated) over-harrier hopping versus quantum- mechanical 
urider-harrier immeVmg. In the classical regime, the rate 
decreases with temperature until reaching the tunneling 
regime, where the corresponding rate approaches a con- 
stant, temperature-independent value. It is worth stress- 
ing that, apart from diffusion problems (in the bulk or 
on surfaces), this typical behavior can also be found in 
many chemical reactions or processes 

According to recent He spin-echo (HeSE) measure- 
ments of H and D diffusion on Pt(lll) [^] at a low cov- 
erage [9 = 0.1 ML) and their ulterior analysis, there is 
a significant contribution of quantum effects for temper- 
atures ranging from 80 K to 250 K, well above the (as- 
sumed) crossover temperature to deep tunneling, Tc ~ 
65 K. As it is shown, quantum motion takes place mainly 
through incoherent tunneling between nearest-neighbor 
sites, which furnishes strong support for a dissipative 
transition-state theory (TST) description. The analysis 
of surface diffusion processes by means of this theory al- 
lows to determine hopping rates, diffusion constants, 
friction coefficients, etc. Therefore, the mobility of ad- 
particles along the surface as a function of the tempera- 
ture and/or the coverage can be closely described. In the 
particular case of Ref. |6|, the authors consider that the 
expression for the hopping rate is given by [loj 



where wq is the vibration frequency at the adsorption 
site, the factor 3 reflects the number of possible jump 
directions, Ea is the barrier activation energy, fc^ is 
Boltzmann's constant, and T is the surface tempera- 
ture. In the Ohmic regime, the expression 2xbU)}j = 
—7 + 4a;^ -|- 7^ gives the friction-dependent barrier fre- 
quency, where 7 is the friction, Xf, is the dimensionless 
Kramers' factor, and ujb gives the vibration frequency of 
the inverted potential at the transition state. On the 
other hand, S accounts for the quantum, temperature- 
dependent correction prefactor [J [l^] , 
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where f„ = 2Trn/h/3 are the so-called Matsubara frequen- 
cies. When applied to the experimental data [6], this 
model renders ujq = 31 mcV for H, LJb/i^a = 5.4 and 
Ea ^ 83 meV for both H and D diffusion, Xb = 0.21 
for H, and xt — 0.29 for D. Because both cob and ojq 
scale as to~^/^, the values corresponding to D can be 
readily obtained by considering the multiplying factor, 
•\/ niH /mo- With these parameters, the crossover tem- 
perature Tc — huJbXb/^TrkB [l] gives Tc ~ 65 K for H 
and Tc « 63 K for D, which indicates that the range of 
temperatures in the experiment is above the crossover 
temperature. According to this model, the conclusion 
arising from the fittings reported in Q is that the hop- 
ping for both H and D is dominated by tunneling, with 
the corresponding rates displaying a certain curvature in 
the region of lower temperatures. However, when the 
friction coefficient 7 is evaluated for each adparticle, it is 
found that 7h ~ 1.15 ns~^ and 70 ~ 0.57 ns^^, which 
display the correct factor 2 between them, but are un- 
expectedly large compared to systems involving heavier 
particles |11|. 

Here we analyze these experimental findings by means 
of an alternative approach based on the so-called bounce 
technique, which allows to determine rates whenever the 
range of parameters (surface temperature and friction) 



2 



are such that the diffusing adparticle tunnels coherently 
and incoherently [l^ 13 1. In particular, in the in- 
coherent tunneling regime, a master equation scheme, 
the so-called Chudley-Elliott model from surface diffu- 
sion [13, [i3> rules the corresponding dynamics. When 
considering the dynamics on a surface, adparticles are 
moving in a two-dimensional periodic lattice of binding 
sites. In the case of activated diffusion, where the ther- 
mal energy is higher than the barrier height, adparticles 
are mainly assumed to perform discrete jumps between 
neighboring sites. Within the Chudley-Elliott model, this 
jump diffusion dynamics on a two-dimensional Bravais 
lattice is accounted for by a master equation in terms of 
the van Hove G'(R, t)-function (or time-dependent pair 
correlation function) [3], widely used to describe a sta- 
tistical ensemble of interacting particles. Thus, given a 
particle at the origin at some arbitrary initial time i = 0, 
G(R, t) gives the averaged probability of finding the same 
or another particle at the surface position R at time t. 
Notice that this function generalizes the well-known pair 
distribution function g(R) from statistical mechanics [13] 
by providing information about the interacting particle 
dynamics. The master equation corresponding to inco- 
herent tunneling can be expressed in general as 



G(R,i) = E^j G(R + j,i), 



(3) 



where Vj accounts for the tunneling rate involved in the 
transition between the lattice point R and the nearby 
one R -|- j, with j being the jump vector among different 
lattice points running over all lattice vectors (positive, 
negative, and even zero). In the case of light particles, 
such as H and D, and for thermal energies lower than 
the barrier height, a tunneling-mediated diffusion pro- 
cess is expected to dominate this nearest neighbors dy- 
namics [6|. 

The space Fourier transform of the G-function is the 
intermediate scattering function, which reads as 



/(AK,i) = (e-'■^K.R(^)g^AK.R(o)^^ 



(4) 



with the brackets denoting an ensemble average. This 
function measures the time correlation loss for a given 
parallel (along the surface) momentum transfer of the 
probe particle, AK. Therefore, it can also provide us 
with information about friction coefficients at different 
coverages and (surface) temperatures along with the ob- 
served AK direction. Without loss of generality, along 
this direction the diffusion process can be assumed to be 
one-dimensional among the different wells formed by the 
surface corrugation, so that the intermediate scattering 
function can then be expressed as a Fourier series as 

/(Ai^||,t)=^G.(t)e^^^ii", (5) 

n 

where Ai^n is a dimensionless momentum transfer re- 
sulting from the projection of the lattice vector j along 



the direction pointed by AK multiplied by the lattice 
constant a, i.e., Ai^n = a||AK|| cos a, with a being the 
angle between AK and j. Here, only first neighbors are 
considered and therefore, given the geometry of the two- 
dimensional lattice for the Ft (111) surface |l|, for four 
of these neighbors lAifnj = aA/'ircos(7r/6), while for the 
other two Ai^n = 0. Regarding n, it labels the n-th well 
of the binding site (bearing in mind this tight-binding 
like model). Now, for nearest neighbors, Eq. ([3]) can be 
expressed in terms of the G„(t) coefficients as 

Gnit) = l^+_iGn-l{t) + V-+^Gn+l{t) - {v+ + i^,7)G„(i), 

, , (6) 

which is precisely the same expression provided by Weiss 
and Grabert for incoherent tunneling in a periodic po- 
tential Il3|. Within this context, the rates v. 



in the 



above expression give the tunneling rate from the (n=p 1)- 
th well to the n-th one, while ± account for the tunneling 
rates to the right or left neighboring well, respectively . 
Tunneling rates are assumed to be equal for the left or 
right direction and independent of the well site. The dif- 
ferential equation (|6]) can be solved analytically for the 
initial conditions, G„(0) — (5„,o- From it, the intermedi- 
ate scattering function is written as 



/(Ai^||,t) = e 



-2rtsin^(Ai<'||/2) _ ^-Vt cos(AKi/^) 



-Vt 



J2 /n(n)e*^^ii". 



(7) 
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where /„ is the modified Bessel function of integer order n 
and T = 1^/2 describes the global tunneling rate. Equa- 
tion d?]) is precisely the function recorded in the HeSE 
technique [6]. 

The bounce technique [§, [H, [13] enables an analytical 
expression for F, which reads as 



F = 



^A^ 

2 uq 



hujr, 



F(C+l/2)- 



(8) 



In this expression, F(-) denotes the Gamma function, A 
is the dressed tunnel matrix expressed as a function of 
the so-called bare tunnel matrix Aq [^ , and C represents 
a dimensionless friction coefficient, given by 



c 



ma 
2'Kh 



1- 



(9) 



Within this analytical model, the diffusion coefficient can 
also be obtained through the relation D = a^f where one 
clearly sees that the usual Einstein relation does not hold, 
since it does not scale with the friction coefficient as 7~^. 

Experimentally, the tunneling (or hopping) rates F are 
obtained after carrying out a best-fit procedure through 
Eq. ([8]) for a given momentum transfer, temperature and 
coverage. In Fig. [U it is displayed the fitting of the ex- 
perimental tunneling rates for H (experiment: blue open 
circles; fitted values: blue solid line) and D (experiment: 
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FIG. 1: Tunneling rates as a function of temperature for H 
(experiment: blue open circles; fitted values: blue solid line) 
and D (experiment: red open triangles; fitted values: red solid 
line) for a coverage 6' = 0.1 ML. 



red open triangles; fitted values: red solid line). In both 
cases, the coverage is 9 = 0.1 ML and temperatures range 
from 80 K to 250 K. Although the quality of the fitting 
is as good as in Ref. here only two parameters were 
needed: a multiplying prefactor, C, and the friction co- 
efficient 7] [through Eq. i.e., 



fit 
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After substituting the corresponding values into Eq. ([9]) , 
we find that the friction coefficients are 7 = 13.6 ps^^ 
for H and 7 = 7.4 ps~^ for D, which are much smaller 
than those previous reported Q and closer to those for 
heavier absorbed species [11|. The prefactors obtained 
are C = exp(32.21) for H and C = exp(32.66) for D. 
From Eq. we also find that the friction coefficient for 
D is roughly a factor 2 smaller than for H, in agreement 
with the mass scaling factor expected for this coefficient, 
■m~^ (isotopic effect). Finally, due to the fact C > 1 
for both H and D (Ch — 2.63 and Cd = 2.85) we would 
like to stress the suitability and validity of this model 
to describe the range of temperatures considered here; 
even valid down to zero surface temperature as stated by 
Grabert and Weiss [l^ . 

In order to analyze the dependence on the surface cov- 
erage, an alternative fitting has also been considered by 
using a temperature-dependent collisional-friction model, 
the so-called two-bath model 18-20|. Within this new 



model, one bath describes the effect of surface phonons, 
while the other one accounts for the collisions among 
the interacting adsorbates. These baths are assumed to 



FIG. 2: Tunneling rates as a function of temperature for H 
(blue solid line) and D (red solid line) for a coverage of 9 = 
0.01 ML. For comparison, the same curves from Fig. 1 are 
plotted for a coverage of 9 = 0.1 ML but fitted to the so- 
called two-bath model (see text). In the inset, the negligible 
variation of the prefactor r((")/r((:^ -I- 1/2) along the range of 
temperatures studied is also shown (solid lines: 9 — 0.1 ML; 
dashed-dotted hue: 61 ^ ML). 



be uncorrelated and therefore the corresponding frictions 
are additive. Accordingly, the total friction, 77, is now a 
sum of two contributions: the usual substrate friction, 7, 
and a coUisional friction. A, accounting for the collision 
among adsorbates (ry = 7 -I- A). For a simple estima- 
tion of the dependence of A on the coverage, one can 
consider Il9l the elementary kinetic theory of transport 
in gases [17| as well as the Chapman-Enskog theory for 
hard spheres. This procedure renders a simple relation- 
ship between the coUisional friction A and the coverage 9 
at a certain temperature T, 



A 



6p9 ksT 



(11) 



where p is the effective radius of the adparticle. Thus, 
by increasing the coverage, the coUisional friction is also 
increased and therefore the total friction. This allows us 
to extend straightforwardly our expressions for quantum 
diffusion by tunneling to different coverages (up to more 
or less 15% for this model). In order to distinguish in 
the friction coefficient the contributions coming from the 
surface and from the collisions, the fitting was carried out 
expressing the exponent in ([8|) as 1 — 67 — c •\/r/1000. 

Following the fitting procedure mentioned above with 
the available experimental data, the curves for 9 =^ 
0.1 ML have again been computed, finding C = 
exp(32.54), b = 0.388 ps, 7 = 13.56 ps-\ and c = 
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0.462 K-i/2 for H, and C = exp(32.93), h = 0.775 ps, 
7 = 7.34 ps-\ and c = 0.382 K^^/^ for D. The newly fit- 
ted curves for 6* = 0.1 ML, displayed in Fig. [5] (blue and 
red dashed lines for H and D, respectively), again show a 
fairly good agreement with both the experiment and the 
fitting curves displayed before (see Fig.[T]). Furthermore, 
the values of the substrate friction 7 are similar with 
those found before assuming that all friction effects were 
included in 77 fulfilling the isotopic effect. Now, accord- 
ing to this new functional dependence of the tunneling 
rate on the coverage and surface temperature, if we con- 
sider c — c'0/0.1, then the curves for a lower coverage 
can also be readily determined, e.g., for Q = 0.01 ML, 
as seen in Fig. [2] (blue and red solid lines for H and D, 
respectively). Moreover, this fitting model also allows us 
to justify our former preassumption that the multiplying 
factor r((^)/r(C-|-l/2) can be considered to be a constant. 
As seen in the inset of Fig. [2l the variation of this func- 
tion along the range of temperatures considered when C, 
is expressed as a function of temperature increases only 
about 1% for H and 0.8% for D as T decreases in the 
case of B = 0.1 ML (for « this increment is mean- 
ingless). This prediction thus shows that a decrease in 
the total friction (e.g., by lowering the coverage or equiv- 
alently the coUisional friction) leads to an enhancement 
of tunneling which is consistent with our model. This 
behavior is no longer valid for higher coverages (greater 
than 15 %) as reported by previous experimental results 
of the same systems by using the helium atom scatter- 
ing technique and analyzed in terms of an Arrhenius law 
[if. Therefore, the application of this two-bath model 
could be extended to higher coverages letting A free for 
the fitting procedure. 

To conclude, the characteristic temperature power law 



I2I . Il3l | for tunneling 



face temperatures, where over-barrier hopping should be 
predominant for heavier adsorbates, one may conclude 
that quantum effects are relevant in the diffusion of light 
adparticles. Apparently, the TST framework does not 
follow such a power law leading to a limited application 
at low surface temperatures. 
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([8]) proposed by Grabert and Weiss 

rates is apparently valid at all temperatures (even zero 
temperature) provided C > 1, as it is in our case. Indeed, 
the application of the Chudley-Elliott master equation 
to describe the incoherent tunneling together with the 
use of such a rate have rendered values for the friction 
parameters in the order of the ps~^, which are within 
the expected range. In general, the loss mechanism may 
come from the lattice relaxation and/or the electronic 
contribution. As reported by Sundell and Wahnstrom 
, this typical power law was first suggested by Kondo 
21j and seems to be related to the electronic contribu- 
tion, i.e., the non-adiabatic response of the conduction 
electrons to the adparticle (H or D) motion. Neverthe- 
less, in view of the temperature-dependent behavior dis- 
cussed elsewhere l^ll3| and also presented here, the use 
of this power law seems to be more general, and valid for 
any loss mechanism or a sum of several uncorrelated loss 
mechanisms (sum of friction coefficients due to electrons, 
phonons, and adsorbates). If the fitting of the experi- 
mental results is quite good in a very wide range of sur- 
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